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Abstract. For certain types of quantum graphs we show that the random matrix 
form factor can be recovered to at least third order in the scaled time r from periodic- 
orbit theory. We consider the contributions from pairs of periodic orbits represented 
by diagrams with up to two self-intersections connected by up to four arcs and explain 
why all other diagrams are expected to give higher-order corrections only. 

For a large family of graphs with ergodic classical dynamics the diagrams that 
exist in the absence of time-reversal symmetry sum to zero. The mechanism for this 
cancellation is rather general which suggests that it also applies at higher orders in 
the expansion. This expectation is in full agreement with the fact that in this case 
the linear-T contribution, the diagonal approximation, already reproduces the random 
matrix form factor for r < 1. 

For systems with time-reversal symmetry there are more diagrams which contribute 
at third order. We sum these contributions for quantum graphs with uniformly 
hyperbolic dynamics, obtaining +2t^, in agreement with random matrix theory. As in 
the previous calculation of the leading-order correction to the diagonal approximation 
we find that the third order contribution can be attributed to exceptional orbits 
representing the intersection of diagram classes. 
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1. Introduction 

The recent work of Sieber and Richter |2] has renewed the hope that spectral 
correlations in systems with chaotic classical analogue can be explained within periodic- 
orbit theory. The universality of these correlations, known as the BGS conjecture [Hj, 
is supported by overwhelming numerical evidence j^- On the other hand there is no 
satisfactory theory for individual chaotic systems, i. e. without any disorder averages. 
Numerically it was found that on time scales longer than the ergodic time of the classical 
analogue, the fluctuations in the energy spectrum of a quantum system follow those of 
an appropriate ensemble of random matrices. For random matrices, the form factor 
K{t), which is the Fourier transform of the spectral two-point correlator, is 

KcoEir) = 2t-t log(l + 2r) (0 < r < 1) 

= 2r - 2t^ + 2t^ + O(r^) , (1) 

for systems with time-reversal symmetry (TR), or 

i^GUE(r) = r (0 < r < 1) (2) 

for systems without time- reversal symmetry (NTR) 

The semiclassical limit of the form factor in a quantum chaotic system can be 
written in terms of a double sum over periodic orbits (PO) using the Gutzwiller trace- 
formula [6.- On short times the relatively small number of contributing periodic orbits 
allows explicit calculation, however the number of POs proliferates exponentially with 
time, so evaluating the sum exactly quickly becomes impossible. In any case the 
universality of the BGS conjecture suggests that beyond the ergodic time the form factor 
does not depend on the specific dynamics of the given system. Berry [7j explained that 
this universality arises from the combined contributions of the huge number of ergodic 
POs. He then calculated the form factor, neglecting all correlations between POs other 
than exact symmetries. Within this "diagonal approximation" , he obtained the leading 
order in r of the random matrix theory (RMT) result. Efforts to reproduce Eqs. (p^, Q 
beyond the diagonal approximation have been limited in their success. At present there 
is no way to derive the series expansion of Eq. from the POs of any chaotic system, 
nor is there a good explanation of why Eq. (0) happens to be exactly reproduced by the 
diagonal approximation for r < 1. 

Currently we only know how to go beyond the diagonal approximation in a 
few special cases. In ^ |2] it was shown, that for uniformly hyperbolic and time- 
reversal invariant billiards on surfaces with constant negative curvature the second- 
order contribution — 2r^ is related to correlations within pairs of orbits differing in the 
orientation of one of the two loops resulting from a self- intersection of the orbit. We 
went on to derive the same result for a large family of quantum graphs jHJin] with ergodic 
classical dynamics, in particular our result was not restricted to uniformly hyperbolic 
dynamics jTHj. A subsequent study [TTj indicated that the mechanism generating the 
contribution — 2r^ also works for systems with antiunitary symmetries other than simple 
time-reversal. 
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Given these results it is a plausible conjecture that in analogy to disordered 
systems ^2] the terms in the power series expansion of K{t) can be identified with 
the contributions of orbit pairs generated by more and more self-intersections. In the 
present paper we will explore this idea for a particular model system: Extending our 
recent paper PO] we will calculate the form factor up to order for a particular family 
quantum graphs. 

This article is organized as follows: In Section 2 we define our model and explain 
how the form factor can be expressed as a double sum over periodic orbits. In Section 3 
this sum is rewritten in terms of diagrams, representing all orbits with a given number 
and topology of self-intersections. Diagrams resulting in a contribution of order are 
considered explicitly. In Section 4 we show that those diagrams which do not require 
time-reversal invariance cancel each other. The summation over the additional diagrams 
in graphs with time-reversal invariance is performed in Section 5, unfortunately here our 
results are limited to a family of graphs with uniformly-hyperbolic classical dynamics. 
Finally in section 6 we explain how we selected the diagrams which give r^- contributions 
by establishing a heuristic rule which predicts the order of a given diagram's contribution 
without an explicit calculation. 

2. Quantum graphs and periodic-orbit theory 

We consider graphs with vertices connected by a total of B directed bonds. A bond 
leading from vertex m to vertex I is denoted by {m,l). For graphs with time-reversal 
invariance it is necessary that for any bond (m, /) there exists also the reversed bond 
(/,m). We do not rule out the possibility of loops, i.e. bonds of the form (m,m). 

The discrete quantum dynamics on a graph are defined in terms of a B x B 
unitary time-evolution operator 5*^^^ which has matrix elements Sl^,], describing the 
transition amplitudes from the directed bond (m, /) to {l',m')\. The topology of the 
underlying graph is refiected in the quantum dynamics because the amplitudes are non- 
zero only if the two bonds are connected at a vertex, 1 = 1'. We choose 

with (T^/^ denoting the vertex- scattering matrix oX vertex /. An explicit example of such 
a graph will be given in Sectional here we keep the discussion as general as possible. The 
phases 0mi are random variables distributed uniformly in [0, 27r] and define for fixed B 
an ensemble of matrices S**-^^ which can be used for averaging. It is possible to interpret 
this ensemble as an infinite energy average for a given quantum graph with rationally 
independent bond lengths jHj. For a unitary operator like S^^^ the form factor is defined 
at integer times t = 0, 1, ... by 

ir(^)(r)=5-i(|tr5f){<^}, (4) 

where r is the scaled time r = tjB. See ||4j for more details on the description of 
two-point correlations for unitary operators. For finite 5, the form factor Q should 

X We drop the parentheses when a bond is used as an index of a matrix. 
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be compared to ensembles of unitary random matrices of dimension B (CUE for NTR, 
COE for TR). However, we are interested here in the hmit of large graphs B —>■ oo, 
keeping the scaled time r fixed 

K{t) = lim K'^^\t), (5) 

B^co 

because this is equivalent to the semiclassical limit of chaotic systems |H]. It is in this 
limit that the form factor is expected to assume the corresponding universal form 
or (0). 

Associated with the unitary matrix S is the doubly stochastic matrix M with 

It defines a Markov chain on the graph which represents the classical analogue of our 
quantum system [HJ Ej- The matrix M can be considered as the Frobenius- Perron 
operator of the discrete classical dynamics. Matrix elements of powers of this operator 
give the classical probability to get from bond (m, /) to bond {k, n) in t steps 



pit) 



M* . (7) 



Under very general conditions it can be shown that the dynamics generated by M is 
ergodic and mixing§, i. e. for fixed B and t — > cxd all transition probabilities become 
equal 

Piii)^ik,n) - as t ^ OO V(m, /), (A;, n) . (8) 

However, since in Eq. (0) the limits S — > oo and t — oo are connected by fixing r, 
Eq. dHJ is not sufficient for showing agreement between PO expansion and RMT. We 
need a stronger condition such as 

^ as 5 ^ oo V(m, /), {k, n) . (9) 

This was already discussed in |15j in connection with the diagonal approximation. In 
fact the precise condition may in principle depend on the order to which agreement 
with RMT is required. In jTOj we derived a condition sufficient for the leading-order 
correction to the diagonal approximation which was slightly stronger than Eq. (jH)): 
The speed of convergence to equidistribution with increasing B cannot be arbitrarily 
slow. However, exponential convergence (corresponding to a spectral gap of M which 
is bounded away from zero uniformly in B) is sufficient in any case. We will restrict 
ourselves to graphs which obey this condition rather than derive a more precise condition 
for the applicability of Eq. Q to the summation of third-order diagrams. 

A connection between the quantum form factor Eq. and the classical dynamics 
given by Eq. (jH)) can be established by representing the form factor as a sum over 
(classical) POs. We expand the matrix powers of S in Eq. (0)) and obtain sums over 
products of matrix elements Sp^p^^p-^p^ ■ ■ ■ Sp^p^^p^p^Sp^p^^p^p^^. Obviously each such product 
can be specified by a sequence of t vertices. Vertex sequences which are identical up 

§ It is plausible to assume that these conditions are satisfied if the underlying graph is connected and 
one excludes special cases like bipartite graphs 
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to a cyclic shift give identical contributions and will be combined into the contribution 
of a periodic orbit P = [pi, . . . ,pt]. For most POs there are t different cyclic shifts. 
Exceptions to this rule are possible if a PO is a repetition of a shorter orbit, but the 
fraction of such orbits decreases exponentially in t. Moreover, if we assume the existence 
of the limit ©, we can approach it through sequences of prime t, which totally excludes 
repetitions. We obtain 

tr5* = t^Ape''^^ (10) 
p 

with Ap = 11^=1 c"p^^i,p^_i and 0p = Y.l=i'Ppi+i,p, (vertex indices are taken modulo t). 
Substituting this into Eq. (jH) we obtain a double sum over periodic orbits 

K^''\r) = ^^(j:ApA*Qe'(f--f^A . (11) 

We can now perform the average over the phases (pmi associated with the directed bonds. 
If the system does not have time-reversal symmetry, all bond phases can be varied 

independently. The total phase of an orbit, 0p, can be written as linear combinations 

tp\ /p\ 

of the bond phases, 0p = J2im''^im yim, where nij^' counts visits of the orbit P to bond 
(m, /). Then we can average over (pi^ using 

'^^iinl^Ulm-n^l^Ul^A =S^^p^^^Q^. (12) 

/ rh, Im Im 

' 'rim 

Thus averaging over all bond phases, {0}, amounts to picking out only those pairs of 
orbits which visit the same set of bonds the same number of times. Therefore the form 
factor for a quantum graph with no time-reversal symmetry (NTR) is 



.(B) . > i 



j(Q) 

, iiiL Im 

Im 



(13) 



Time-reversal symmetry implies symmetric vertex-scattering matrices 

rr^') = a^^^ (U) 

m'm mm' \ J 

and that the phases of a pair of bonds (m, /) and {l,m) related by time-reversal obey 
(pmi = (t>im- The average in Eq. (fTTjl runs over all independent phases and results in 

l>m 

Thus averaging over all independent bond phases, {(f)}, amounts in this case to picking 
out only those pairs of orbits which visit the same set of bonds, or their time-reverses, the 
same number of times. Hence, the form factor for a quantum graph with time-reversal 
symmetry is 



^2 r 

J_J L 7^ '"^ 



(16) 

,^ Im ' ml ' Zm ' ml 

' l>m 

If the graph is defined to have bonds with fixed lengths and magnetic vector potential as 
in jH] , we can average over an infinite energy window (for NTR-systems we also average 
over the vector potential). Then the orbit pairs contributing to the form factor are those 
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where P, Q have exactly the same lengths. On a graph with rationally independent bond 
lengths this is equivalent to 

NTR: nl^^=nlS yi,m (17) 

TR: nj^) + = + n!?? V/,m, (18) 
so this averaging procedure also leads to Eqs. (fT^ and (fT^ . 



3. The expansion in self-intersections of the periodic orbits 
3.1. From orbits to diagrams 

The calculation of the form factor is now reduced to a combinatorial problem: The sum 
over the pairs P, Q in Eq. (fT^ and Eq. (fT^ must be organized such that Eqs. (fTTj) or 
(|T8|l are satisfied. This can be done by composing P and Q from the same segments, or 
arcs, which appear in P and Q in different order and/or orientation. This is possible if 
the orbit P contains self-intersections, i.e. vertices which are traversed more than once, 
see Fig. ^ for examples. In general an orbit P has many self-intersections and many 
partner orbits Q satisfying Eqs. (fTTj) . (fTHj) such that a summation over all possible Q for 
a fixed P is too complicated. Instead we fix a permutation of arcs followed by the time 
reversal of selected arcs and sum first over all possible pairs of orbits P, Q related by this 
transformation. The clearest way to represent all possible transformations is graphical 
(Fig. Q), hence we refer to them as diagrams. The sum over all diagrams finally gives 
the form factor. 

The main problem with this approach is to ensure that each orbit pair P, Q is 
counted once and only once. This is difficult because for some pairs P, Q the operation 
transforming P into Q is not unique. Such orbit pairs are relatively rare in number but 
nevertheless they give essential contributions to the form factor lOj. We will explain 
our techniques for preventing the double counting of orbit pairs in Sections l3 .21 and 

If we consider P as a single arc with no intersections, Q = P is the only possibility 
in the NTR case. For TR the orientation of the arc can be reversed such that Q = P is 
a second option. The corresponding diagrams have a simple circular shape. Summation 
over these orbit pairs is nothing other than the diagonal approximation. It produces 
-f^NTRi = T and -R'tri = 2r, respectively. In \Ui^ we considered orbits, P, made from 
two arcs, 1 and 2, joined at a single intersection a and evaluated the (off-diagonal) 
contribution corresponding to the resulting 8-shaped diagram. We found this gave rise 
to the second-order term in Eq. (^, Ktr2 = — 2r^, while there is no contribution of this 
order for a NTR-system. In this paper we calculate the r^-contribution by assuming 
that P contains three or four arcs connected at intersections. A discussion of why only 
these particular diagrams contribute to the r^-contribution is deferred to Section IHl 

We begin by listing in Fig. ^ all diagrams which contribute at third order in r. We 
denote arcs by numbers 1,2,... and the intersection points by Greek letters a, P, . . .. An 
arc can be identified by a sequence of vertices, which does not include the intersection 
vertices, or, alternatively, by a sequence of bonds, which includes the bonds from and to 
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NTR3a 

orbits : lQ'2/?3a4/3 and la4/33a2/3 
factors : 1/4 

restrictions : (1 — 5.5, s^) (1 — Ss^,,^) 



NTR3b 

orbits: Ia2a3a and laSa2a 
factors : 1/3 

restrictions : (1 — Ss^a-,) (1 ~ ^ais-,) ( 



TR3a 

orbits : lQ'2/33a4/3 and la4/32a3/3 
factors : 1/2 

restrictions : (1 - 63234) (1 - ^s^fi) 
with 1 7^ 2 and 1 7^ 3 

TRSb 

orbits : la2l3?,(3Aa and lai(}Zj32a 
factors : 1/2 

restrictions : (1 — 5.,, 74) (1 — ^,54/,) 
with 1 7^ 1 and 3 7^ 3 





Figure 1. Topology of NTRSa, NTR3b, TR3a, TRSb and TR3c. In each case a pair 
of orbits is shown. One follows the solid line throughout (in the direction marked by 
triangular arrows). The second follows the solid lines (possibly with reversed direction) 
except at the intersections (denoted by circles) where it follows the dotted line. Each 
circle represents a single vertex where a self-intersection of the orbit occurs. Next to 
each topology we give the corresponding weight factor ('Section l3.2|l and the restrictions 
f Section 13. 3|) . Restrictions are indicated by the double-headed arrows. Solid arrows 
indicate a "full" restriction of the form (1 — (5^ ), while dotted arrows indicate a "half" 
restriction of the form (l — ^(Jy). 
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the intersection points. The length of the ith arc is denoted by tj and is defined as the 
number of bonds in the arc (which is one more than the number of vertices in the arc). 
The sum of the lengths of all arcs gives t, the length of the orbit. The length of an arc 
is at least one. For an arc i leading from a to /3 we denote the first vertex following a 
by Si and the last vertex before /3 by /«. In the degenerate case when the arc going from 
a to /9 is the single bond (a,/?) and does not contain any vertices (tj = 1) this implies 
Si = P and fi = a. 

As shown in Fig. ^ the arcs forming an orbit P and its partner Q are identical, 
but the way they are connected at the intersections differs. The orbit P is given by the 
connections drawn as continuous lines, while its partner orbit Q is given by connections 
drawn as dotted lines. The orbits P and Q are also written as a symbolic code to the 
left of each diagram: a path that goes from the beginning of arc 1 to vertex a then 
on arc 2 to vertex f3 and so on is denoted as Ia2f3 ■ ■ ■. The diagrams in Fig. ^ divide 
into two classes, NTR and TR. In the NTR-diagrams all the arcs of Q have the same 
orientation as the corresponding arcs in P, while in the TR-diagrams some of the arcs 
of Q are time-reversed. For a system with no time-reversal symmetry, only the two 
NTR-diagrams are possible, thus there are two r^-contributions to the form factor, 

-^NTRS = -K^NTRSa + -^NTR3b • (19) 

For a system with time-reversal symmetry (TR), diagrams in both classes contribute 
and the form factor is a sum of five terms 

KtR3 = 2 (-ft'NTRSa + -f^'NTRSb + -f^TR3a + -f^TR3b + -f^TRSc) • (20) 

The factor of two is due to the fact that for every diagram in Fig. ^ there is another one 
with Q replaced by its complete time- reversal, Q, which gives an identical contribution. 

3.2. Avoiding double- counting I: multiplicity factors 

The set of diagrams possesses certain degeneracies which can be accounted for by simple 
prefactors multiplying the contributions. One such degeneracy is taken care of by 
the factor of two in Eq. (j20|l . In this subsection we discuss how to determine the 
other multiplicity factors arising due to the cyclicity of the POs and symmetries in the 
diagrams. To sum over all orbit pairs P, Q for a given diagram we sum over all possible 
arcs forming the orbit P. Consider the diagram NTR3a as an example. Let li and I3 
be some fixed arcs starting at (3 and ending at a, while I2 and ^4 denote arcs from a to 
(3. As we sum over all possible realisations of arcs 1,2,3,4 in NTR3a, we encounter a 
particular orbit P where these arcs are given by 

l = /i, 2 = k, 3 = /3, 4 = ^4. (21) 

However we also encounter the orbit P' where the arcs are 

l = /3, 2 = /4, 3 = /i, A = k. (22) 

The orbit P' is related to P by a cyclic shift and, therefore, it is actually the same 
orbit. As we are focusing on pairs of orbits, we check the partner orbits resulting 
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from P and P', too. The partners for P and P' are Q = [/i, a, /4, /3, /a, a, ^2, /3] and 
Q' = [l3,a,l2, PJiyCtJi, P], respectively, and they are also related by a cychc shift. 
Hence in the process of summation we will encounter the same orbit pair four times, 
once for each of the four possible cyclic permutations of P. To compensate for this we 
introduce a multiplicity factor of a quarter. 

To put this formally, we denote by q{P) the operation transforming P into Q for 
a given diagram, e. g., gNTR3a([1234]) = [1432] and gTR3a([1234]) = [1423]. Further we 
denote by a the (cyclic) left shift of the symbolic code, i. e. (t([1234]) = [2341]. To 
determine the multiplicity factor we need to count all cyclic permutations a'' such that 

qoa^P) = a''' oq{P) or g o (t'=(P) = a'^' o ^(P) (23) 

for some k' and arbitrary P, i. e. application of q to the shifted code yields the same or 
the completely time-reversed result, up to another shift. Obviously the second option 
is only applicable to the TR diagrams. Noting that the trivial solution k = k' = is 
always available, we proceed to list the factors for each diagram. 

• NTR3a: We have gNTRSa ocr'' = a^~^ o gNTRSa for any /c = 0, 1, 2, 3 and consequently 

^NTRSa = 4. 

• NTR3b: Similarly we have gNTRSb o = a^^^ o gNXRSb for k = 0,1,2 and 

^NTRSb = 3. 

• TR3a: The only nontrivial solution to Eq. pHj) is = 2 where gTRSa o = gxRaa- 
Therefore we have mxRSa = 2. 

• TR3b: The only nontrivial solution is k = 2 where gxRSbO o"^ = cr^oQ'TRSb- Therefore 
we have mxRsb = 2. 

• TR3c: Eq. has no solution besides the identity k = 0, therefore mTRSc = 1- 

When we evaluate the contribution from each diagram we will include a factor 1/m in 
order to compensate for the ambiguity just described. 

3.3. Avoiding double- counting II : restrictions and exceptions 

As shown in Ref. ^01; tangential self-intersections of orbits are a potential source for 
double-counting of orbits which must carefully be avoided. By a tangential intersection 
we mean the situation where an orbit does not merely cross itself but follows itself for 
at least one bond. For example the orbit 

> fi a 13 Si+i > fj+i 13 Sj+i ^ ■■■ (24) 

crosses itself along the non-directed bond {a, (3). It is easy to mistakenly count such an 
orbit once with a as the intersection point and once with (3 as the intersection point. 
We avoid this using a method outlined in Ref. JUI- We uniquely define the intersection 
point by ruling that if there is an ambiguity then the intersection is as far to one side 
as possible. As an example we show some ambiguous intersections in Fig. |21 and insist 
the intersection is as far to the right as possible. For the 2-intersection we do this by 
demanding that vi 7^ f2, this is achieved by introducing a factor of the type (1 — S^^-o^), 
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Figure 2. Examples of ambiguous intersections, where we have removed the 
ambiguity by placing the intersection (the shaded vertex) as far to the right as possible. 
This is enforced in the left-hand diagram by introducing a factor of (1 — (Jtjj^j), while 
in the other diagrams it is enforced by a factor of (1 — Sy-^y^Sy-^y^). 

referred to as a restriction on the diagram. For the 3-intersections a restriction of the 
form (1 — Svj^y^Sy-^v-i) removes the ambiguity. However, we will not actually use this 
restriction on any 3-intersection, because we will see that stronger restrictions apply in 
the diagrams we evaluate. 

For NTR3a, TR3a and TR3b we choose the following restrictions to ensure the 
ambiguities at intersections are removed 

• NTR3a: Antrsr = (1 - ^^2^4) (1 " ^^1^3) 

• TR3a: AxHSa = (1 - ^sisj (1 - ^sau) 

• TR3b: ATR3b = (1 - 5s,u) (1 - 

where Sj denotes the first vertex on the arc i and fi denotes the last. We wish to 
emphasize that there is no unique way of imposing the restrictions, since they are merely 
convenient ways of excluding the double counting of certain contributions. What is more, 
the individual results for NTR3a, NTR3b, TR3a, TR3b and TR3c may depend on the 
particular choice of restrictions. Only the final sums in Eq. ()19p and ((201) ^'^^ depend 
on them. 

Now that we come to NTR3b, we will see exactly how much freedom we have in 
choosing restrictions. First we want to ensure that we count tangential intersections 
correctly. For a 3-intersection, such as the one in NTR3b, we could do this by 
setting ANTRSb = (1 ~ ^siS2^s2S3)- However we also notice that there are ambiguous 
contributions which could be counted in either NTR3a or NTR3b |: 

• NTR3b with either si = S2, S2 = S3 or si = S3 is equivalent to NTR3a with any of 
the following: (t^ = 1 & = /g), {t^ = 1 k f2 = f^), (ta = 1 & /i = /a) or {t^ = 1 

& /2 = /4) 

• NTR3b with either /i = /2, /2 = /a or fi = f^ is equivalent to NTR3a with any of 
the following: (ti = 1 & Si = S3), (^2 = 1 & S2 = S4), (t^ = 1 Sz Si = S3) or (t^ = 1 

& S2 = S4). 

Obviously we should only count each of these contributions once, but we have the 
freedom to choose whether we count each of them in NTR3a or NTR3b. The physical 

I This abundance of choice, when any of the three NTR3b diagrams is equivalent to any of the four 
NTR3a diagram is another manifestation of the cyclic symmetry discussed in Section 13.21 and is taken 
care of by the multiplicity factors. 
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quantities (|19p and (j2(Jj) contain the sum of NTR3a and NTR3b, so all choices are strictly 
equivalent. However given that we have imposed the restriction si 7^ S3 on NTR3a the 
second type of orbits (NTR3b with fi = fj) cannot belong to NTR3a. Thus, once we 
have chosen the above restrictions for NTR3a we are forced into the choice 

ANTRSb = (1 - f^sisa) (1 " ^sisg) (1 " ^3283) ■ (25) 

Before we can move on to TR3c, we must first look carefully at the restriction 
we placed on TR3a. In Section f3.2l we introduced the factor of 1/2 to avoid double- 
counting. The double-counting in this particular instance was caused by the permutation 
cr^ = [3412] which swaps around arcs 1^3 and 2 4 and produces a pair P' = o"^(P) 
and Q' = cr^(Q), which is identical to P, Q up to a shift. However the restriction S2 7^ f^ 
that we introduced on TR3a is not symmetric with respect to this permutation. For 
the orbits satisfying si 7^ /2 and S3 7^ f^ this does not present any problems. Let us 
consider what happens when arcs 1 and 2 are different but have si = f2- This orbit is 
still counted twice in the summation over all possible arcs, but in the second instance 
the intersection point /5 is shifted, resulting in t[ = t^ + 1, t'2 = t^ + 1, t'^ = ti — 1, 
t'^ = t2 — 1- We illustrate that by the following example of orbits which contribute to 
TR3a. The pair 

P = [P,'j,a,a,b,'j, P,d,a,c] and Q = 7, a, a, c, /3, 7, 6, a, (i] (26) 

is obtained by combining the arcs 

arc 1 = {(3, 7) — > (7, a) (a, a) arc 3 = {(3, d) — > (ci, a) (27) 

arc 2 = (a, 6) (6, 7) ^ (7, /?) arc 4 = (a, c) ^ (c, /3) (28) 

with the intersection points a and /?, or by combining the arcs 

arc 1 = (7, P) — i> (/?, d) — > ((i, a) arc 3 = (7, a) — ^ (a, a) (29) 

arc 2 = (a, c) (c, /3) -> (/3, 7) arc 4 = (a, 6) ^ (6, 7) (30) 

with the intersection points a and 7. We therefore see that the factor of 1/2 works also 
when Si = f2- But not when, in addition to si = /2, ti or t2 is equal to 1. These orbits 
appear in the sum for TR3a only once and are subsequently multiplied by 1/2, so it 
appears we miscount their contribution. On the other hand these orbits can also be 
counted in TR3c, as shown below. We find it convenient to keep the above restriction 
on TR3a, thus counting half their contribution in TR3a. This forces us to count the 
other half of their contribution in TR3c. 

Now we can move on to finding the restrictions on TR3c. First we list the special 
cases of TR3c which could be counted in other diagrams. 

(i) TR3c with si = f2 is equiv. to TR3a with (t^ = 1 Sz S3 = f^) or (ti = 1 & si = /2) 

(ii) TR3c with S2 = fs is equiv. to TR3a with (^4 = 1 & S3 = /4) or (t2 = 1 Sz si = /2) 

(iii) TR3c with S2 = S3 is equiv. to TR3a with (ti = 1 & /i = f^) or (^3 = 1 & /i = f^,) 

(iv) TR3c with fi = f2 is equiv. to TR3a with (t2 = 1 & S2 = S4) or (^4 = 1 & S2 = S4) 

(v) TR3c with si = /i is equiv. to TR3b with (t2 = 1 & ^2 = fi) or (t4 = 1 & S4 = /2) 
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(vi) TR3c with S3 = /a is equiv. to TR3b with (^2 = 1 & S2 = fi) or (^4 = 1 & S4 = 

Now we carefully count in TR3c only those contributions which have not already been 
counted in TR3a or TR3b. Lines (ji)) and (jUl) above show that cases si = /2 and S2 = /s 
should be counted in TR3c with the factor 1/2. Line (jm)) shows that the case S2 = S3 
should not be counted in TR3c as it is fully counted in TR3a; the case /i = /2, line 
(jiv|l . should be fully counted in TR3c. Lines (jvj) and show that the cases si = /i 
and S3 = /3 are not counted in TR3b and should be fully counted in TR3c. All this is 
realised by the restrictions 



Above we have ensured that no orbits are double-counted among the diagrams 
NTR3a, NTR3b, TR3a, TR3b and TR3c. However we should also exclude the orbits 
that have already been counted at lower orders of the expansion. Considering NTR3a, 
if arc 1 is identical to arc 3 (or arc 2 identical to arc 4), the diagram is reduced to 
giving a contribution to the diagonal approximation, so it should not be counted here. 
Fortunately, the restrictions we have put on NTR3a ensure that this contribution is not 
counted. Moving on to NTR3b, if any two arcs in the NTR3b diagram are self-retracing 
the diagram reduces to a diagram already counted as a r^-contribution in a TR-system. 
Therefore, in the TR case, we should subtract its contribution from the sum. However, 
we will see at the end of Section Wl\ that such a contribution is zero. 

For TR3a, we insist that 17^2, 17^3, 47^2 and 4 7^ 3 because the orbits breaking 
these rules have already been counted at C[r^] of the expansion. For the same reason 
we insist that TR3b obeys 17^1,37^3, 27^4, while TR3c obeys 1 7^ 1, 2 7^ 2 and 
3 7^ 3. Note that some of the restrictions are superfluous since they refer to orbits that 
are already excluded. For example we can drop the restriction 4 7^ 3 because this is 
automatically enforced by the stronger restriction S3 7^ /4. 

The complete set of restrictions is the following 



• NTR3a: AwTRSa = (1 - Ss.sJ (1 - Ss,s,) 

• NTR3b: ANTRSb = (1 - S^.s,) (1 - Ss,s,) (1 - Ss,s,) where orbits with (2, 3) = (2, 3) 



• TR3a: ATR3a = (1 - Ss,s,) (1 - ^^3/4) with 1 ^ 2 and 1 ^ 3. 

• TR3b: ATR3b = (1 - 5.2/4) (1 - 5.4/2) with 1 ^ 1 and 3 ^ 3. 

• TR3c: Atr3c = (1 - 5.,sJ (l - l^s^/a) (l - ^5,,f,) with 1 ^ 1, 2 ^ 2 and 3 7^ 3. 



We reiterate that this self-consistent set of restrictions is not unique. And, although 
this choice lead to simpler calculations than all the others we tried, we can not rule out 
the possibility that there is another self-consistent set of restriction which would further 
simplify our calculations. 

3.4- Orbit amplitudes 

Before we can attempt the summation over all orbit pairs P, Q within a given diagram, 
we still need to understand the structure of the product ApAn appearing in Eqs. (fT!?jl . 




(31) 



must be subtracted for systems with TR symmetry. 
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We consider the diagram NTR3b as an example. Let arc 1 be of length ti, 
consisting of the vertices [xi,X2, ■ ■ ■ Xt^^^i], where Xi = si and Xt^^^i = fi. Then both 
Ap and Aq will contain factors cr^.^^^, "^igfii; ^xT,x2^ ■ ■ ■ ' '^(^^^^i-l ■ Thus when we evaluate 
the product ApAg, the contribution of the arc 1 will come in the form 

\'^X2,a'^X3,xi'^X4,X2 <Ja,xt^-2\ — ^{a,xx)^{xx,X2)^...^{xt^~i,a) — -t^l , 1<J^J 

which is the c/a55^ca/ probability of following arc 1 from bond (a, si) to bond (/i,a)§. 
Analogous considerations for arcs 2 and 3 lead to the probabilities P2 and P3. The 
factors not yet accounted for in Pi, P2, P3 are the transition amplitudes picked up at the 
intersection vertex a: 

ApA*^ = p,xP2xp,x aiy:\ai:i X {ai:},aiy:iy. (33) 

To evaluate the contribution of a given diagram a product like Eq. (jHHj) must be summed 
over all free parameters, namely all intersection points and all possible arcs connecting 
these points. The latter summation includes a sum over the lengths ti of these arcs with 
the restriction that the total length of the orbit is t. 

The summation over all the intermediate vertices 0:2, ^3, . . . ,Xt^-2 along arc 1 can 
be performed immediately, because it is unaffected by the restrictions discussed in the 
previous subsection. This summation adds the classical probabilities of all possible paths 
leading from bond (a, Si) to bond (/i, a) in ti — 1 steps and results consequently in the 
classical transition probabihty P(asiy-^(fi a) given by Eq. Analogous summations 
over the other arcs produce P'(^^ and P'^^^ ^^^^(^f^ The above approach extends 
trivially to the TR diagrams when we recall that time-reversal symmetry implies that 
the matrices o"^^^ are symmetric. 

The remaining summation is over the lengths tj of all arcs, the first and the last 
vertex and fi of all arcs i with ti > 1 and the intersection points like a. For 
general graphs this sum is still too complicated for explicit calculations, mainly because 
transition probabilities like P^^^'^^^^^j^ are dependent on the details of the topology 
of the graph. For sufficiently long arcs, however, these transition probabilities can be 
replaced by B^^ according to Eq. Q. Then the sum over vertices decouples into a 
product of sums associated with each self-intersection vertex a which can finally be 
evaluated using the unitarity of the vertex-scattering matrices 0"*^°^ . This is the strategy 
we shall follow in the next two sections, where explicit summation of the NTR3 and 
TR3 diagrams is performed. 

4. Summing the NTR contributions 

4-1. Summation of NTR3 diagrams 

Starting with the NTR3a diagram, we write 
1 1^ 

^NTR3a(^) = J-^^^\t - Ei=l ij] ^NTR3a X PNTRSa X AwTRSa , (34) 

§ Pi = 1 if arc 1 contains no vertices, i. e. if ti = 1. 
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where 



P _ p(*l-l) p(*2-l) p(t3-l) p(t4-l) /0(^\ 

^NTRSa - ^(^,,,)_(/,,„)^(„,s2)->{/2,/3)M/3,^3)^{/3,a)M">S4)^(/4,/3) 

ANTRSa = (1 - 5siS3)(l - '^S2S4)- (37) 



of generahty we assume that ti > t/4. From Eq. we have ^-j ~ -B ^ and 

the only factors in Eq. (fH^ depending on /i are crl"/! "^Ir/i ■ Using the unitarity of the 
(T-matrices we perform the summation 

E-;i-s:= 5^2.4- (38) 
/i 

However the restriction ANTRSa contains the term (1 — 6s2S4,), leading to the result 

i^NTRSa = . (39) 

Calculation of -ft'NTRSb goes along the same route with 

^NTR3b(r) = ^ J2 J2 SNTR3b X PNTR3b X ANTR3b , (40) 

where 

2.NTR3b - f^.3/2 ^^2/1^^1/3^^2/3 ^^3/1 ^si/2 (41) 

P _ p{*l-l) p(*2-l) p{*3-l) {Ar^\ 

^NTR3b - ^(Q,si)^(/i,a)^(a,S2)^(/2,a)M",S3)^(/3,a) ^^"^1 

ANTR3b = (1 - 5siS2)(l - <^S2S3)(1 - 5s3sJ. (43) 

Exactly as for NT3a we can sum over where arc i is long, which results in 5-function 
which we combine with ANTR3b to get the answer 

i^NTR3b = . (44) 

The sum of the NTR3a and NTR3b diagrams vanishes and so 

i^NTR(r)=0. (45) 

Thus we see that for a wide class of quantum graphs without time-reversal symmetry 
the r^-contribution to the form factor is zero, as expected from the BGS conjecture. 

We wish to note that the derivation given above is relatively simple, since NTR3a 
and NTR3b both vanish due to the particular choice of restrictions which make orbit 
pairs in the intersection of NTR3a and NTR3b unique (see the discussion near Eq. (j25j) ). 
Had we chosen to assign all ambiguous diagrams to NTR3a, then the results for NTR3a 
and NTR3b would both have been non-zero, although the total sum K^'j^^^t) would of 
course still have equalled zero. 

To apply the above result to the TR case, described by (pUj) . we must subtract 
the contribution of the NTR3b diagram with two self-retracing circs, cLS discussed in 
Section Eini We use the fact that long self-retracing arcs give only exponentially small 
corrections because the number of free summation variables is reduced by a factor of 
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Figure 3. The picture on the left is a self-intersection in a NTR-contribution with 
seven crossing arcs. If two bonds leaving this intersection coincide, in this case 2 and 
6, the intersection can be redrawn (as on the right) as more than one intersections. In 
this case there are three intersections, one 2-intersection, one 3-intersection and one 
4-intersection, the latter two being at the same vertex. 



two. Thus we only need to consider short self-retracing arcs. Without loss of generality 
we can assume that arcs 2 and 3 are self-retracing and short, implying that ti must be 
long enough for p(*i-i) = B-^ to hold. Then the sum over fi results in a factor of 6s2S3 
and, combining this with the restriction ANTRsb, we find that the contribution of this 
case is identically zero. Thus the NTR diagrams contribute nothing to the form factor 
of TR-systems. 

4- 2. Generalization to higher orders 

One may speculate that the arguments given in the previous subsection admit a 
straightforward generalization to higher-order diagrams. Given an nth-order diagram, 
we impose the following restriction on each of its intersections 

^ = I[{^-Ks,) , (46) 

where the product is over the set of all arcs leaving the intersection. Now we can 
evaluate diagrams in the same way as we did for n = 3. As soon as <^ t/terg at least 
one arc must be long. If arc i is long then p(*'~i) ~ B^^ and the sum over fi generates a 

5- function. Combining this (5-function with the restriction at the vertex produces zero. 

To justify the choice of the restriction ()4(j|l for any intersection, we notice that if 
any two bonds leaving the vertex are the same, the intersection can be rearranged as a 
group of more than one intersection, each satisfying the above restriction. An example 
of such rearrangement is presented in Fig. IHl If the original intersection was part of an 
nth-order diagram, then the rearranged one is part of another valid nth-order diagram 
(as can be shown by counting the powers of B, see Section IH)). The above restriction 
thereby helps to prevent double counting of orbits with tangential intersections. 

This argument essentially shows that the contribution of all nth-order diagrams is 
zero in the NTR case. However, an important detail is missing: one has to show that 
all eligible pairs of periodic orbits belong to one and only one diagram, i.e. that we 
did not miss anything and did not count anything more than once. Unfortunately we 
found pairs of orbits that violate both parts of this statement. These counterexamples 
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Figure 4. A fully-connected graph with four vertices and sixteen directed bonds. 



seem to be "rare", in the sense that the sum of their contributions vanish as — oo, 
however a rigourous proof of this observation remains an open problem. 

To summarise, a generahzation of the argument of Section 14.11 sketches a proof 
of exactness of the diagonal approximation for r < 1 in the absence of time-reversal 
symmetry. To complete the proof one would have to verify that the above restriction 
counts all relevant pairs of orbits once and only once. 



5. Summation of TR3 diagrams for a fully-connected "Fourier" graph 

Evaluating i^TRs for a general class of graphs is a complicated and tedious task P^ . 
Fortunately, the calculation simplifies considerably for a special case described below. 
In this section we restrict our attention to fully-connected graphs with N vertices and 
B = N"^ directed bonds, including a loop at each of the vertices. An example with 
= 4 is shown in Fig. 0] We assume that the vertex-scattering matrices are 

(I) 1 / mm'\ , , 

= exp 2vn^ (47) 



for all /. These matrices were proposed in J3] and result in a particularly fast 
convergence to RMT-like statistics. Because of the analogy to the discrete Fourier 
transformation from m to m', we call a vertex endowed with the scattering matrix ()47j) 
a "Fourier" -vertex. The corresponding matrix M represents uniform scattering at the 
vertex, 

M^>^^ = \a^2'rn? = , (48) 
and thus the probability to get from (a, b) to (c, d) in t step is 



(a,b)^(c,d) V J{d,c),ib,a) \ B~\ t>l 

(a,6)^(b,a)' 



It is also useful to have an expression for P(ab)^(ba)^ probability to get from (a, b) to 



{b, a) following only self-retracing paths. The contribution of each path is iV~* and, due 
to their special structure, a self-retracing path of 2m + 1 or 2m -|- 2 steps will contain m 
vertices not including the initial a and b. Each of these m vertices can be freely chosen 
from the vertices of the graph, resulting in N"^ different self-retracing paths. Thus, 
the probability -P^^^^)^^^ takes the form 

P{a,b)^{b,a) = N-^~", with t = 2m or 2m + 1 , (50) 
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i. e. it is indeed decaying exponentially in time. 
5.1. Summation of TR3a 

Here we calculate the contributions of orbits with the topology of TR3a which obey 
the conditions S2 7^ S4 and S3 7^ f^. We enforce these conditions by multiplying the 
contribution of all orbits of this topology by 

A^n3. = {l-Ss,s,){l-Ss,u). (51) 
Thus the contribution of TR3a is 

i^TR3a(r) = ~J2^[^- ^t=M E E ^TRSa X PxRSa X AxRSa, (52) 



^ ^ {U} a,f3 Siji 



where 



_ ^(") ^(") (/3)* W* 

P _ p(*l-l) p(*2-l) p(*3-l) p(t4-l) /t-/|\ 

^TRSa - ^(/3,s,)_(/i,o)^(„,^2)_^(/2,/3)^(/3,^3)^(/3,c,)^{a,S4)^(/4,/3) " l^^^^l 

If arc 1 is long enough to be ergodic, the sum over si simplifies to 

(1 - ^s,u) E ^S4<i: = (1 - '^.3/4)^.3/4 = . (55) 

SI 

If arc 2 or arc 3 are ergodic we can carry out a similar sum over /2 or f^ respectively, 
these sums also yield the answer zero. Thus we can only get a non-zero contribution 
when arc 4 is the only ergodic path. However, we argue in Section IHl that we need at 
least two arcs to be long (ergodic) since otherwise the contribution can be neglected. 
Thus the non-zero contribution discussed above will only give a small correction which 
will vanish in the limit B oo. The two restrictions 1 7^ 2 and 1 7^ 3 do not change 
the above argument at all, so we have ignored them. We conclude that 

A'TR3a(r)=0. (56) 

5.2. Summation of TR3b 

Here we calculate the contributions of orbits with the topology of TR3b which obey the 
conditions S2 7^ /4, S4 7^ /2, 1 7^ 1 and 3 7^ 3. The first two conditions will be enforced 
by means of a factor 

ATR3b= (l-5s2/4)(l- 5^4/2), (57) 
the latter two we will enforce below "by hand" . Thus 

A'TR3b(r) = ~J2^[^- ^t=M E E STR3b X PTR3b X AxRSb , (58) 



^ {U} a,f3 SiJi 



where 



^TR3b - f^S4/3^.3/2^S2/l^.l/4 ^/2 f3^^3i4^/4/l^^li2 V^^) 



's4/3"^3/2"s2/l^Sl/4 hh"" S3Si^ hh'' S1S2 
P _ p(*l-l) p(*2-l) p(*3-l) p{*4-l) 

^TR3b - ^(„,s,)^(/,,o)^(c,,s2)^(/2,/3)M/3,^3)^(/3,/3)M/3,S4)-(/4,a) " ^^"^^ 
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We only need to consider cases where ti > 3 and ts > 3 because shorter arcs are purely 
self-retracing (1 = 1) and so must be excluded. We will treat the restrictions 1 7^ 1 
and 3 7^ 3 using the following inclusion-exclusion procedure: the sum in (|HH|l with these 
restrictions is equal to the sum without the restrictions, minus the sum with 1 = 1, 
minus the sum with 3 = 3, plus the sum with both 1 = 1 and 3 = 3. 

The first sum yields zero after the summation over si or over S3 in a fashion similar 
to Eq. ()55|) . The second sum we perform with respect to S3 while the third is summed 
with respect to Si, in both cases the answer is zero. Thus K^RshiT) is equal to the sum 
with both 1 = 1 and 3 = 3, which can be written as 

{ti} 

X E E ^£s;H{.i,a)^S!s;)l(s3,/3) X ^TR3b X ^TRSb X ^siji^ssja, (61) 

where we used the fact that P'^^^^^^i^^^ = B^^ for j = 2, 4, while -P(^^^^)^(^ is defined 
above eq. (jSHl)- Upon substitutions /i = Si and f^ = S3, and using the symmetry of the 
matrices a, Eq. (fT^ . STR3b becomes 

T — |2|^(o) |2|^(/3) |2|-.(/3) |2 _ Ar-4 /^.qN 

We also notice that the probabilities P do not depend on the start and end bonds. Now 
we can perform the sum over a, (3 and all Sj and /j, which, taking into account various 
delta-functions, gives the factor N^{N — 1)^. We get 

i^TR3b(r) = I'^^^^^^^f^ E ^ - P''^-''P''^-'\ (63) 

^ ^ {U} 

from which it is clear that the dominant contribution comes from ti = 3, 4 and ^3 = 3, 4; 
the contributions from other values of ti and ^3 are of order 0{N~^). After carrying out 
the sum over ^2 using the 5-function which forces t^ = t — t2 — n with n = ti + t^ = 6,7,8 
we get 

1 ^2 t-3-n .3 

i^TR3b=4x-- E 1 = V = 2^' ^^^^ 

where we have dropped corrections which vanish in the limit B, N ^ 00 and the factor 
4 comes from the number of possible choices of ti and ^3. 



5.3. Summation of TR3c 

Here we calculate the contributions of orbits with the topology of TR3c which obey the 
following restrictions. First we should only count half the contribution when S2 = fs 
or si = /2. Secondly S2 7^ S3, 1 7^ 1, 2 7^ 2 and 3 7^ 3. The restrictions which apply 
to whole arcs will again be enforced "by hand" using an inclusion-exclusion procedure 
similar to that used above, the rest of the restrictions are 

AtR3c = (1 - Ss^ss) (1 - \^S2h) (1 - l-^si/a) ■ (65) 
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KTRs^ir) = ^ E - ^i=iU] E E Strsc X Ptrsc x Atrsc, (66) 



{U} O. SiJi 



where 



V — ^(") rr(") 

^TR3c - (^s^h^^s^h'^s.h ^/2/3^^3/i^sii2 . l67j 

P _ p(*l-l) p(*2-l) p(t3-l) ff^Q\ 

-TtRSc - -r(a,si)^(/i,a)-^(a,s2)^(/2,a)Ma,S3)^(/3,«) " 

The summation here is similar to the sums in TR3b: first we ignore the restriction 1 7^ 1 
(but enforce the restrictions 2 7^ 2 and 3 7^ 3) and carry out the sum over /i to get 

(1 - ^...3) E ^S/.^S: = (1 - Ss,s,)Ss,s, = . (69) 
/i 

Then we subtract the sum over orbits with 1 = 1 (again enforcing the restrictions 2 7^ 2 
and 3 7^ 3). Similarly to TR3b, it turns out that the dominant contribution comes 
from orbits with ti = 3,4, i.e. p(*i~^) = N~'^. Since ti <^ t we use the argument from 
Section El to note that we are only interested in orbits where ^2,^3 ~ t and thus both 
arc 2 and 3 are ergodic. This leaves us with 

i^TR3c(r) = -2 X E E E StR3c X AtR3c X 5,, j„ (70) 

where n = 3, 4. We sum over /i using Eq. (fT^ to get 



^TR3c 



I ^ 2 

S2Sl 



'^Si/3'^S3/2^/2/3^*3^1 - ^Sl/3 ^53/2^/2/3^5331' l^'^^l 

open up the brackets in Atr3c, 

^1 - \^S2f'^ (1 - I'^^i/a) = (1 - i<^Sl/2) - i'^S2/3 + i^52/3'^5l/2) (^2) 



and are now facing the sum 

E ^ll^ll^S^S: (1 - ^S2S,) [(l - i5../2) - i5s2/3 + ¥s2f.A.f.] ■ (73) 

a,Si,S2,/2,S3,/3 

Invoking the unitarity of the cr-matrices, it is an easy exercise to show that this sum 
evaluates to N^{N - l)/2 + N{N - l)/4. 

Combining the above information and ignoring subdominant contributions we arrive 

at 

i^TR3c(r) = -2-^E^'(^- l)/2 = -r^- (74) 

<2 



5.4. The TR3 result 

Remembering that we proved (NTR3a + NTR3b)=0 in Section ITTl we simply need to 
substitute the results of the three previous subsections into Eq. pO|l to get 

i^TR3(r) = 2r=^ . (75) 

Combining this result with the one in ^0] proves that the form factor for the fully- 
connected Fourier graph coincides with the GOE form factor up to the third order in 
r. 
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6. Estimating the order of a diagram 

In this section we discuss a rule for finding all diagrams which contribute to the nth 
order in the small r expansion of the form factor. The rule is 

(#of arcs) — (#of intersections) = {n — 1). (76) 

Thus for n = 2, we need only one diagram which is (2, 1) in the format (#of arcs, T^of 
intersections), and this is the contribution we considered in ^U]. Here we are interested 
in n = 3, so we must consider both (3, 1) and (4, 2). It is these diagrams that we show 
in Fig. HJ 

To get the rule ()76p we count powers of i? in a diagram's contribution. Equations 
()13|) and |T?)|l have a pref actor of so a r"-contribution to the form factor must 
get from the summation over the orbits. In the ergodic limit, according to 

Eq. 0, each arc will contribute the weight B~^, while each intersection contributes the 
weight B, thus we have Eq. (fTBjl . The origin of the factor of B associated with each 
intersection can be explained as follows. First of all, the set of all vertices {vj} adjacent 
to an intersection point 7 can be split into two equal subsets satisfying the following 
property: if there is a transition (wj,7) (7, I'm) in either P or Q then Vj and Vm 
belong to different subsets^. This is particularly simple for the NTR3 diagrams where 
the two sets are simply {sj} and {/«}• If we now do the summation over all vertices in 
one subset and invoke the unitarity of the scattering matrix at the vertex 7, the result 
will be a product of 5-functions SuiU2^u2U3 ■ ■ -^Ukui where Uk are the vertices from the 
second subset, ordered in an appropriate way. Now the summation over U2, ■ ■ ■ ,Uk will 
give 1 while the summation over ui and 7 will give the sought-after factor of B, since 
the only restriction on Ui and 7 is that they have to be two ends of the same bond. 

To make this recipe work for the diagonal term (~ r^), the corresponding diagram 
being just a looping arc, we need one extra ingredient, the starting vertex for the loop. 
The position of the vertex is not determined, it can be placed anywhere on the looping 
arc, unlike the intersection points in other diagrams. To compensate for this ambiguity 
when we sum over all periodic orbits fitting such diagram, we divide the sum by the 
number of vertices in the loop, t. 

Now we discuss why counting of powers of t does not work for obtaining Eq. ()76|) . 
Let us estimate the order of t for a given diagram. Firstly, there is in the prefactor 
of Eq. (fT^ or (fT^ . Secondly, for a diagram with a arcs, the lengths tj of arcs satisfy 
J2iti = t thus the sum over all possible ti gives a factor proportional to Then of 
the diagrams in Fig. ^ NTR3b and TR3c appear to have four powers of t while the 
rest have five. Similarly, the diagram we evaluated in ^U] gets three powers of t. The 
leading contributions to all the diagrams appear to have at least one more power of t 

I In other words the graph buUt on vertices Vj, connected if there is a transition {vj,j) — > {-j,v„i), is 
bipartite. This graph is nothing else but the structure drawn inside the circles in Fig. ^ The graph is 
bipartite since it is 2-regular (the valency of each vertex is 2) and each connected component contains 
an even number of bonds. 
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than they should§. However, we show in ^U] and the present article that the numerical 
coefficient of this "out of order" term is zero — at least for diagrams contributing up to 
third order in r. 

The arguments given above in favour of Eq. (f7B|) are certainly too vague to be 
considered a proof. In particular, we cannot presently check our assumption that 
terms giving incorrectly large powers in t disappear also for more complicated diagrams. 
However, summing rigorously the contributions of all the diagrams obtained from this 
set of rules we show a posteriori that we indeed get an expansion which depends only 
on the scaled time r. We also take confidence in our method from the fact that our rule 
generates the same diagrams that were used in perturbative calculations of the form 
factor for disordered systems with the non-linear sigma model |12j . 

Nevertheless, counting powers of t is very useful in the following situation: If for 
some reason the lengths of some arcs are forced to be fixed, the estimated power of 
t can drop low enough that we can safely ignore the contribution of such a diagram 
in the B —>■ oo limit without actually evaluating it. In particular, we see that to get 
a no n- vanishing contribution of order t^, at least two arcs in any diagram must have 
unrestricted lengths. Note that this does not mean that there is no contribution from 
orbit pairs where the maximum length of all arcs is restricted. For any given B and t 
there are orbit pairs with so many self-intersections that the maximum arc length is less 
than the time required for ergodicity. Then the method discussed in this paper must 
fail, this may explain why the power series expansion in r breaks down at r = 1 for 
NTR (and at r = 1/2 for TR) despite the fact that the PO sums Eqs. (jlHj) . (jlfij) are 
exact. 

7. Conclusions 

At first sight, the achievements of the present paper may appear moderate. What is the 
point in going from second to third order in a series which is infinite, in particular if this 
step becomes possible only by restricting the range of systems considered to very special 
models? But we believe that such a point of view is too short-sighted. Semiclassical 
theories are obviously indispensable for a complete understanding of spectral statistics 
in systems with chaotic classical dynamics. As they are not based on a random matrix 
conjecture, such theories have the potential to account for important system-specific 
corrections which one can hope to extract once the emergence of universality within 
semiclassics is fully understood. On the other hand, the restriction to the diagonal 
approximation has so far severely limited the success of the semiclassical approach. 

Going beyond the diagonal approximation in semiclassical PO theories is possible, 
as demonstrated by Sieber and Richter PJ Ej. But when more than a ffist-order 
correction is required, one will inevitably encounter the problems discussed in this paper. 
For example, one needs methods to select diagrams contributing at a given order, and 

§ This, however, is not the case for the "diagonal" diagram where the power of t is right, which explains 
why advancing beyond the diagonal approximation was (and still is) so hard. 
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we suggested a solution in Section IHl It will be necessary to account for the ambiguity 
introduced by the representation of the form factor in terms of diagrams, and we have 
solved this problem at least in quantum graphs for the diagrams which contribute up 
to third order fSection l3.3p . It is important to note that a variety of orbits give the 
relevant contributions, and our calculations in Section indicate that the generalization 
of the leading-order correction to higher orders cannot be achieved by considering a 
single type of diagram only. 

Here the third-order result for TR systems is limited to a class of uniformly 
hyperbolic quantum graphs. However, we have no reason to believe that any conclusion 
will be substantially changed when the calculation is done for a more generic model as, 
e. g., in [T^. While going further than third order for TR-systems is beyond us at the 
moment, the prospects for doing this in NTR-systems are more promising. We hope 
that the method presented in Section will prove applicable there. 

8. Acknowledgments 

We gratefully acknowledge extremely productive discussions held with M. Sieber and 
U. Smilansky. We also would like to thank the group of Prof F. Haake in Universitat 
Essen for hospitality and exchange of ideas. We are grateful to a referee who found 
a large number number of misprints and inconsistencies in the paper and helped us 
improve the presentation of Section lT^ A significant part of the research was performed 
while both GB and RW were working at the Weizmann Institute of Science, Rehovot, 
Israel. During that time GB was supported by the Israel Science Foundation, a Minerva 
grant, and the Minerva Center for Nonlinear Physics; and RW was supported by the 
U.S.-Israel Binational Science Foundation (BSF) and the German-Israel Foundation 
(GIF). HS thanks the Weizmann Institute for kind hospitality during a number of visits. 
RW is now supported by an EPSRC Grant. 

References 

[1] M. Sieber and K. Richter. Physica Scripta, T90:128, 2001. 

[2] M. Sieber. J. Phys. A, 35:613, 2002 

[3] O. Bohigas, M. J. Giannoni, and C. Schmit. Phys. Rev. Lett., 52:1-4, 1984. 

[4] F. Haake. Quantum Signatures of Chaos. Berlin: Springer, 2000. 

[5] M. L. Mehta. Random Matrices Boston: Academic Press, 1991. 

[6] M. C. GutzwiUer. J. Math. Phys., 12:343-358, 1971. 

[7] M. V. Berry. Proc. R. Sac. Lond. A, 400:229-251, 1985. 

[8] T. Kottos and U. Smilansky. Phys. Rev. Lett, 79:4794-4797, 1997. 

[9] T. Kottos and U. Smilansky. Ann. Phys., 274:76-124, 1999 

[10] G. Berkolaiko, H. Schanz, and R. S. Whitney. Phys. Rev. Lett., 88:104101, 2002 

[11] P. A. Braun, F. Haake, and S. Heusler. J. Phys. A, 35:1381, 2002 

[12] R. A. Smith, I. V. Lerner, and B.L. Ahshuler. Phys. Rev. B, 58:10343, 1998. 

[13] F. Barra and P. Gaspard. J. Stat. Phys., 101:283-319, 2000 

[14] J. R. Norris. Markov Chains, Cambridge: Cambridge University Press, 1997. 

[15] G. Tanner. J. Phys. A, 34:8485, 2001 



Form factor of quantum graphs 

[16] G. Berkolaiko. Preprint: |nrin.CD/0305009| 2003. 



